Symmetry-Enhanced Performance of Dynamical Decoupling 
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We consider a system with general decoherence and a quadratic dynamical decoupling sequence 
(QDD) for the coherence control ol a qubit coupled to a bath of spins. We investigate the influence 
of the geometry and of the initial conditions of the bath on the performance of the sequence. The 
overall performance is quantified by a distance norm d. It is expected that d scales with t, the 
total duration of the sequence, as r mln { Jv x,Jv i }+i^ wnere ]\j x anc [ j\[ z are the number of pulses of the 
outer and of the inner sequence, respectively. We show both numerically and analytically that the 
state of the bath can boost the performance of QDD under certain conditions: The scaling of QDD 
for a given number of pulses can be enhanced by a factor of 2 if the bath is prepared in a highly 
symmetric state and if the system Hamiltonian is SU(2) invariant. 

PACS numbers: 03.67.Pp, 03.65.Yz, 82.56.Jn, 03.67.Lx 
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I. INTRODUCTION 

Improvements both in resonance spectroscopy and in 
quantum information rely on the ability of suppressing 
unwanted couplings between the system and its envi- 
ronment. Uncontrolled couplings arc often the origin 
of phase accumulation and in general of decoherence. 
Therefore, a faithful manipulation and preservation of 
quantum states is required. 

The dynamical decoupling (DD) is an open-loop con- 
trol scheme to average out the undesircd coupling be- 
tween the system (qubit) and the environment (bath) by 
means of stroboscopic pulsing of the qubit. The DD was 
developed by Viola and Lloyd [l[ from the original idea 
of Hahn 0. 

In its original formulation the DD makes use of equidis- 
tant 7r pulses to average out only a single coupling along 
one spin direction, usually the z direction, (pure dephas- 
ing) - we think for example of the Carr, Purcell, Mai- 
boom and Gill (CPMG) sequence [!, H|. A remarkable 
advance is the optimal DD discovered by Uhrig Q , whose 
sequence has the minimun number of pulses for a given 
order of the suppression of the decoherence. It was shown 
that UDD can also be used to suppress longitudinal re- 
laxation . Recentl y ot her non-equidistant sequences 
have been proposed pl-QjJ . 

The most general case concerns the suppression of de- 
phasing and longitudinal relaxation at the same time. A 
sequence of pulses having a single level of suppression 
cannot suppress general dephasing. Sequences with two 
sorts of pulses have been pro pos ed where concatenated 
sequences are used, like CDD \l2\ and the CUDD Q for 
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example. 

Recently West et al. [l3| have proposed a near opti- 
mal scheme that suppresses arbitrary couplings to order 
t n (t is the duration of the total sequence) between the 
qubit and the bath using 0(N 2 ) pulses. The sequence 
consists of two levels of nested UDD, therefore the name 
quadratic UDD (QDD). The validity of UDD can be ex- 
tended to analytically time-dependent Hamiltonians [T3| 
which is an important ingredient for the demonstration 
of QDD. 

Wang et at showed that the effect of QDD can 
be decomposed in the effects of the inner and the outer 
sequences. The concept of mutually orthogonal opera- 
tion set (MOOS) for nested Uhrig DD was introduced: a 
set of control operators on the inner level is not affected 
by a set of control operators on the outer level if both 
sets come from a MOOS. Higher order protection of a 
MOOS can be achieved if even-order UDD sequence on 
different levels are nested. Thus the results in Ref. [TH 
demonstrate the validity of QDD with even-order UDD 
sequence on the inner level. If the inner level has an odd 
order the symmetry group generated by MOOS is broken 
and the scheme based on nested UDD can not be applied 
anymore. It appears that this problem has been solved 
by Jiang and Imambekov (lfjj who have provided a proof 
of the validity of nested UDD (NUDD) sequences that 
relies on a mapping between NUDD and a discrete quan- 
tum walk in 2m dimensional space. The case of QDD 
corresponds to m = 1. At last, an alternative proof of the 
validity of QDD and a numerical investigation of the scal- 
ing of the errors along specific spin directions for QDD 
has been presented in Ref. [l?] and in Ref. [TH, respectively. 

In this paper we want to draw the attention to the ef- 
fects of the state of the bath on the performance of the 
sequence. The fact that the specifics of the bath can 
limit the performance of a sequence is already known 
[ijj [2(| . It was tested experimentally that UDD can out- 
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perform CPMG if the environment is characterized by a 
hard cutoff [H, [2l[ while for soft cutoffs equidistant se- 
quences perform either better or the same [22j. Other- 
wise UDD seems to perform very well for electron spins 
in irradiated malonic acid crystals (2~l| as well as for ap- 
plications of magnetic resonance imaging [23j . 

So far we have always viewed the environment as an 
unavoidable restraint on the prolongation of the coher- 
ence of a spin (qubit). Hence one has to eliminate or 
at least to reduce the coupling between environment and 
system because the coupling between environment and 
system transfers disorder from the environment to the 
system. But does the environment's disorder always act 
against coherence in the system? 

Here we show that the performance of a given sequence 
can be enhanced if the system Hamiltonian H is SU(2) 
invariant and if the initial state of the bath, i.e., its den- 
sity matrix pe is completely disordered: pe oc 1. We call 
such a state an infinite-temperature state. We simulate 
the effect of a QDD sequence on a bath of spins both for 
a completely anisotropic as well as for an isotropic cen- 
tral Hcisenberg spin model. For both cases we analyze 
the scaling of QDD when the bath is prepared either in 
a product state or in an infinite-temperature state. Four 
cases are studied as summarized in Tab. [TJ The cases, 
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TABLE I: The scaling of QDD with the duration of the se- 
quence t classified according to the degree of symmetry of the 
Hamiltonian and of the initial state of the bath. The notation 
-/Vmin refers to the minimum number of pulses of the inner and 
the outer sequence. 



where the bath state is characterized by a low degree of 
symmetry, independent of the degree of symmetry of i?, 
provide the lower bounds for the scaling of QDD: For 
short times QDD scales always as t"™" 1 " 1 , where N m - m 
is the smallest number of pulses, either of the inner or of 
the outer sequence [13, E3 • In the off-diagonal case in 
Tab. [TJ where H is of low symmetry and the bath state 
of high symmetry, the scaling exponent \ depends on the 
number of pulses. Otherwise, the scaling of QDD is en- 
hanced to the power T 2 ( NsaiXi+1 > if H and pe are highly 
symmetric. 

The paper is set up as follows: in Sect. |H] and in the 
Sect. lIIII the numerical results for the low and for the high 
symmetry cases. In Sect. IIVI we provide the analytical 
argument for the appearance of the factor 2 in the high 
symmetry case. In Sect.|V]we study the off-diagonal cases 
with mixed symmetry. At last we draw our conclusions 
in Sect. ED 



II. CASE 1: LOW SYMMETRY 

We start from the case where both the system Hamil- 
tonian as well as the initial state of the bath have a low 
deg ree of symmetry. We consider a central spin model 
[2J-[28[ characterized by a completely anisotropic Hamil- 
tonian of the form 



H = H 



B "I" -ffqB 

M M M 

= EE^^+E* 05 -^ ( 2b ) 

i— 1 j>i i—1 

where all nine entries of the 3x3 Jq and J\ matrices 
are random numbers drawn from the interval [—1, 1] (sec 
Fig. [TJ). The system does not show any symmetry. The 

entries for the matrices Jq and J{ are fixed randomly 
at the beginning of the simulation and they remain the 
same for all the numerical results we present in this ar- 
ticle. The spin labelled with zero represents the qubit 
while M defines the number of spins in the bath. The 
scaling appears to be essentially independent [29| of M; 
we considered for our simulation M = 8. Calculations 
for M = 3 yield the same results as far as the exponents 
are concerned. 
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FIG. 1: (Color online) Central spin model described by Eq. 
pjl with M = 4. The solid lines represent the coupling be- 
tween the central spin (qubit) and the spins of the bath, while 
the dashed lines represent the couplings between the spins of 
the bath. 

The QDD sequence is made of an outer sequence of 
7r-pulses about a x and of an inner sequence of 7r-pulses 
about a z . The number of pulses for each sequence is N x 
and N z , respectively, and the total number of pulses of 
the sequence is N x + N z + N X N Z . The switching instants 
are given by 



tj = r sin 



2(N X 



2{N Z 



(3a) 



(3b) 
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for j = {!,.., N x } and k = {1, .., N z }. We use the nota- 



tion t 



N x +1 



We start with an initial density matrix of the total 
system of the form 



(4) 



with 7 = {x, y, z}. The first factor in the tensor product 
refers to the Hilbcrt space of the qubit, the second to the 
Hilbcrt space of the bath. Furthermore, we introduce the 
notation 



(5) 



For the low symmetry case we assume that the bath is 
initially in a pure product state 



life) = (X) hi) 



(0) 



so that pb = |V'b)(?/'b|- For the high symmetry case we 
choose 



p B = Ib/D 



(7) 



where D is the dimension of the Hilbert space of the bath 
so that TrePB = 1- 

The overall performance of the sequence is given by 
the norm distance |30| 



d 2 = I E 4 

l£{x,y,z} 

4 := Tr q [A^(r) 



(8a) 
(8b) 



with 



A H (r) := Tr B (u*P - U(r,0)p^U(r,0^ 



(9) 

The norm distance measures the distance of the real evo- 
lution to the ideal one. The operator P is defined by 



a, cr x a z ~<T x ...cr 



(10) 



It incorporates the effects of the pulses. The operator 
U(t, 0) represents the evolution operator of the system 

U(T,0) = T{e- i So*Wt\ (ii) 

where T stands for the time-ordering and 

U B ■■= 1 ® e~' lTliB (12) 

is the dynamics of the isolated bath. During the appli- 
cation of the sequence the sign in front of the coupling 
terms between the qubit and the bath perpendicular to 
the pulse direction changes every time a 7r pulse is applied 



because (XiCTj = —UjUi for i ^ j. Thus, in the toggling 
frame, the system Hamiltonian ([5]) can be written as a 
time dependent Hamiltonian 



m 



H(t) = H B +Y, E ( J l 



(13) 

t=i j,k=i 

where the switching functions fj(t) are a piecewise con- 
stant functions with values ±1. By cry' we refer to the 
k component of the vector (?W . Analogously, we refer by 
(Jl)j k to the (j,k) element of the matrix J\. 



Our simulations show that for N x ^ 0, N z ^ and 



for N x =N z = N^doc t n+1 + O (t n+2 ) 
for N x > N z -)• d cx t n ^ +1 + O (r w =+ 2 ) , 
for N T <N z ^-d(x t n * +1 + O (t n *+ 2 ) . 



(14a) 
(14b) 
(14c) 



The data agree 
for the overall error. For ei- 



as shown in Tab. [TT| and in Fig. 2(a) 
with the results of Ref. 
ther N x = or N z = (UDD sequence) we find that d 
scales as r, as expected for a Hamiltonian with general 
decohcrence. On the other hand, if N z = for example 
and H q B is a pure dephasing Hamiltonian, i.e., no cou- 



plings with a x ^ or with <jy'' occur in the Hamiltonian, 
the norm distance scales as o?udd oc t n " +1 . 

Here we present the results of a given random configu- 
ration of the entries of Jq and J{. We also checked that 
the scaling is the same for other random configurations. 



N Z \N X 





1 


2 


3 


4 


5 


6 





1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1 


1.00 


2.01 


2.00 


2.01 


2.01 


2.01 


2.00 


2 


1.00 


2.01 


2.97 


3.01 


3.01 


3.01 


3.02 


3 


1.00 


2.00 


2.99 


4.01 


4.07 


4.09 


4.09 


4 


1.00 


2.00 


3.00 


4.01 


4.99 


4.98 


5.05 


5 


1.00 


1.95 


3.02 


3.95 


5.01 


5.98 


6.05 


6 


1.00 


1.98 


3.01 


3.94 


5.00 


5.84 


6.95 



TABLE II: Low symmetry case: Scaling exponent £ of the 
norm distance d(r) with r, the total duration of the sequence. 
The reported numbers are determined from the slope of the 
curve d vs. r in a double logarithmic plot. 



Different choices of the initial bath state pb of the 
form ([6]), i.e., varying the |7i), can affect the scaling of 
d Xl d y and of d z . but not the overall scaling of d: If 
d~f oc t" 4 "' then the leading order of the norm distance 
scales as d oc r mm 7{ A T}. Hence the scaling exponent 



reads Q = min 7 {A 7 }. This is what the analytic ar- 
guments require for QDD p^Hl~7| . Hence the analytic 
bounds on the exponents are sharp for the low symmetry 
case. 



III. CASE 2: HIGH SYMMETRY 

We consider an SU(2) invariant isotropic central spin 
model with Hcisenberg couplings. We choose a Hamil- 
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tonian of the form of Eq. <j2j) with Jq = aXj^ 1 and 
J\ = where a and A are two generic constants 

while ]q and j\ are random numbers between —1 and 1. 



N Z \N X 





1 


2 


3 


4 


5 


6 





1.99 


2.00 


2.00 


2.00 


2.00 


2.00 


2.00 


1 


2.00 


3.99 


3.95 


4.00 


3.98 


4.00 


4.00 


2 


2.00 


4.00 


5.99 


6.13 


6.00 


5.99 


5.99 


3 


2.00 


3.99 


5.99 


7.98 


7.97 


8.01 


8.01 


4 


2.00 


4.00 


5.99 


7.99 


9.97 


9.92 


9.94 


5 


2.00 


4.00 


5.99 


7.99 


9.97 


11.93 


11.94 


6 


2.00 


4.00 


6.00 


7.97 


9.95 


12.01 


13.95 



TABLE III: High symmetry case: Scaling exponent C, of the 
norm distance with r for an isotropic Hamiltonian with gen- 
eral decoherence and for an initial density matrix such as in 
Eq. (|15|l . The scaling exponent are derived from a fit of the 
numerical curves for d vs. r in a double logarithmic plot. 



If the bath at t = is described by the following density 
matrix 



Pb oc 1 B , 



(15) 



the suppression of the decoherence is enhanced by a fac- 
tor 2. The simulation for QDD yields the scaling expo- 
nents reported in Tab. lIIII We deduce the following rules: 
for N x ^0 : N z ^0 and 

N x = N z = N -»• d oc r 2 ^ +1 ) + O ( T 2N + 3 ) , (16a) 
N x >N z ^dtx t 2 ^ +1 ) + O (r 2JV *+ 3 ) , (16b) 
N x <N z ^doc r^+D + O (t 2N *+ 3 ) , (16c) 

If N x = or N z = the norm distance scales as r 2 . A 
graphical representation of the data is provided in Fig. 



2(b) 



The state of Eq. (|T5|) is a completely disordered state 
where no particular spin direction or state is singled 
out. Such a state can be referred to as an "infinite- 
temperature" state. This is not unusual in NMR ex- 
periments where already at room temperature one finds 
huj-^jkoF m 10 -5 , where wl is the Larmor frequency of a 
spin and &b is the Boltzmann constant. This means that 
the thermal energy exceeds all internal energy scales by 
many orders of magnitude. 

Note that for ps oc 1b the norm distance J5) 
with the partial Frobenius norm distance I2S 



coincides 



IV. SU(2) INVARIANCE 

The appearance of the factor 2 can be explained in 
terms of the different parity of Hb and H q B under spin 
rotations. We write the Hamiltonian ([2]) in the form 



H = t s <E>A + ' 

li£{x,y,z} 



A,. 



(17) 




(b) 



FIG. 2: (Color online) The data of Tabs. HQ and [HI] are graph- 
ically represented in panel (a) and in panel (b), respectively. 
The numbers are rounded to their first digit. Here we have 
introduced the notation here d a r^, where £ stands for the 
scaling exponents. 



In Eq. (fl~7| the operators Aq and A^ act only on the bath 
while Is and act only on the qubit. Since the identity 
operator and the Pauli matrices form a complete basis 
for all system operators the evolution operator can be 
expanded according to 

U(t,0) = 1 s ®B {t) + (18) 

H=x,y,z 

where Bq(t) and B^ir) are non-trivial functions of the 
operators Aq and A^ and of the switching functions /^(i), 
see Eq. (|A2[) in the Appendix and Ref. |3l|. For the sake 
of simplicity we will omit the time dependence of the 
operators Bq and B^ from now on. From the unitarity 
of U(t, 0) we conclude 



B B t 



E 

fj,e{x,y,z} 



B„Bl 



and 



; E ■ 

fj,,v^{x,y,z} 



(B Bl + h.c. 



(19a) 



(19b) 



for fixed k € {x, y, z} and e M!/K being the Levi-Civita sym- 
bol. In (|19b|) we omitted the non-singular factor propor- 
tional to Pauli matrices because the vanishing must be 
ensured by the bath operators. In the Heisenberg picture 



the density matrix (|1]) evolves according to 



ti\r) = U(r,0)p^U(r^ 

= f/(T,0)(p^ ) ®p B )c/(T,0) t . 



(20) 
(21) 
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We trace out the bath and use the unitarity of U(t, 0) 
([1511 to obtain 



Tr B P £\t)=T^+T^ 



T. 



(7) 



T 



(t) 



with 



T 



(7) 



=4 7) Tr B p B 



T. 



(7) 



E/ (^Ti 

.(7) 



Tl 



(7) 



d^b 



T« := 



E 

-« E 

/i,f,KG{x,y ,2} 



The coefficients 



Tr B [SoPbS/J 



depend only on the bath operators while 



and tf 



(7) 



(7) 



PgCTft are functions of the qubit 
operators only. Note that for a pure dephasing model 
[3lT ] the terms T^ 7 ' and T4 do not appear. 

We consider a global operator P„ that rotates all the 
spins of our system around the v = x, y, or z axis by the 
angle tt. Here wc are interested in the SU(2) invariant 
Hamiltonian such as the one discussed in Sec. IIIII Then 
wc have 



(22) 

(23a) 
(23b) 
(23c) 
(23d) 



(24a) 
(24b) 

(7) 



P V B P} = B P V B V P$ = B V 



and therefore 

6„ = Tr B [BoPbBI] 



= Th 



Tl'T 



P v B Pl P„ PB Pt P v B\Pl 
B P vPB Pl M 



(25a) 
(25b) 

(26a) 
(26b) 
(26c) 



for (i / i/. Thus, if p B is invariant under rotation P„, 
which is the case for p B oc 1b, we can conclude from 
(|26c[) that b^ = for p £ {x, y, z}. In fact, the analogous 
argument also implies b^ .„ = for p ^= v, though we will 
not use this fact here. The condition p ^ v is needed to 
ensure that we can flip the sign of the two factor B^ and 
B\ separately. 

If the coefficients b^ vanish only the terms proportional 
to b^v remain in Eq. (|2"2")l . We know from the analytic 
properties of the QDD sequence that the operators B^ 
with p £ {x,y,z} all scale at least with 7- > min+1 where 
N m [ n :— min(N x , N z ) [l5l - [T7| supported by numerical re- 
sults in Rcf. [l8| and in the present work. Hence the coeffi- 



cients b^ scale with t n " 



:n+1 and the coefficients b^ v with 



T 2Af min +2 jj cnce ^he vanishing of the 6 P terms in Eq. (|22|) 
automatically reduce the decoherence by doubling the ex- 



1+2 with the 



ponent in the scaling d oc t n " 
total duration r of the sequence. 

Note that for a model of pure dephasing, e.g., only 
tJz appears in (|2b[) . wc do not need the symmetry with 
respect to two operators P^. It is sufficient to have either 
P x or P y which invert a z so that we can conclude that 
b z vanishes in order to know that the scaling exponent 
doubles. This was already seen in the numerical data 
presented and analyzed in Ref. [29|. 



CASE 3: MIXED SYMMETRY 



(a) 



N Z \N X 





1 


2 


3 


4 


5 


6 





0.99 


1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


1 


1.00 


2.00 


2.00 


2.00 


2.00 


2.00 


2.00 


2 


1.00 


1.99 


3.01 


3.00 


3.00 


3.00 


3.00 


3 


1.00 


2.00 


3.00 


4.00 


4.00 


4.00 


4.00 


4 


1.00 


2.01 


3.00 


4.02 


5.00 


5.00 


5.01 


5 


1.00 


2.00 


3.00 


4.00 


5.02 


6.00 


6.00 


6 


1.00 


2.00 


3.00 


3.98 


5.00 


6.00 


6.99 



(b) 



N Z \N X 





1 


2 


3 


4 


5 


6 





1.98 


2.00 


2.00 


2.00 


2.00 


2.00 


2.00 


1 


2.00 


2.00 


2.00 


2.00 


2.00 


2.00 


2.00 


2 


2.00 


3.99 


6.00 


5.75 


5.99 


5.99 


5.99 


3 


2.00 


4.00 


5.80 


4.00 


4.00 


4.00 


4.00 


4 


2.00 


4.00 


6.00 


7.96 


9.96 


10.27 


10.30 


5 


2.00 


3.97 


5.97 


7.97 


9.85 


6.00 


6.00 


6 


2.00 


4.00 


5.99 


7.99 


9.95 


11.93 


13.92 



TABLE IV: Mixed cases: In Table (a) the Hamiltonian is 
SU(2) invariant while the bath is initially prepared in a prod- 
uct state. In Table (b) the Hamiltonian is asymmetric of the 
form in Eq. (0, the entries of and J{ are all randomly 
chosen, and pe oc 1b. 



Here we analyze the off-diagonal cases of Tab. [I] They 
are characterized either by a SU(2) invariant Hamiltonian 
and a low-symmetry bath state or by a low-symmetry 
Hamiltonian with a high-symmetry initial bath state 
oc 1b- Because the Hamiltonian and the density ma- 
trix have a different degree of symmetry the analytical 
argument of Eq. (|26a[) for b^ does not hold anymore. 

The numerically found scaling exponents are reported 
in Tabs. IIVI and depicted in Fig. [3l The coupling con- 
stants used to derive the scaling exponents of this Table 
are the same as those used for the simulation of the SU(2) 
invariant Hamiltonian or the ones used for the asymmet- 
ric Hamiltonians used in Sees. [TT1 and Hill 

If the Hamiltonian is SU(2) symmetric and the bath is 
initially prepared in a product state the scaling exponents 
look the same as those of Tab. [TTJ If the Hamiltonian 
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FIG. 3: (Color online) The data of IIVI (b) is graphically rep- 
resented. The numbers are rounded to their first digit. The 
notation £ refers to the scaling exponents of the norm distance 
d oc t . 



is asymmetric and pe oc 1b we find for N x = N z = 
N that the scaling exponent is N + 1 for TV odd and 
2(N + 1) for N even. For N x > N z the exponent is 
either N z + 1 if N z is odd or 2(N Z + 1) if N z is even 
while for N x < N z we find that d scales as 2(N X + I). 
We cannot provide any explanations for this alternating 
behavior of the scaling exponents and for the reason why, 
for N x > N z , it depends only on the number of pulses 
of the inner sequence. This is still an open question for 
future investigation. 

Furthermore, it is worth mentioning that peculiar 
choices of the couplings imply non-generic behavior. If 
the couplings j\ is the same for all i the QDD sequence 
behaves like a UDD sequence with N z pulses for N x odd, 
as if the qubit were subject only to pure dephasing: Ei- 
ther A x ^ or A y ^ in Eq. [17] It is not clear why 
the QDD behaves like a UDD sequence in this case. A 
possible explanation is that the product state we take as 
initial state of the bath is not the most general one. It is 
not entangled. In the Appendix we analyze the first three 
cumulants of the evolution operator for the case N x = 1 . 
We find that they contain only one qubit operator, cither 
a x or <jy, and that the dephasing term proportional to 
a z is always zero. For N x even we recover the results of 
Table HI 



VI. CONCLUSIONS 

We investigated the influence of the state of the bath in 
suppressing general decoherence by means of a QDD se- 
quence. The performance of the sequence is measured by 
the norm distance d which is essentially the norm of the 
remaining decoherence. Thus the performance is quanti- 
fied by the scaling of d with r, the total duration of the 
sequence. Recent papers fla - tlq proved the properties of 
QDD and clarified the dependence of the scaling on the 
number of pulses N x and N z : The overall scaling of QDD 
is given by £ := m\n{N Xl N z } + 1 independent of the de- 
tails of the environment, i.e., d oc t*. In this sense QDD 
is a universal sequence for general decoherence such as 
UDD is a universal sequence for pure dephasing. 

In the present work, we have shown that the actual 



performance of QDD can be even better than expected 
on the basis of the general mathematical arguments. This 
improvement occurs if the Hamiltonian and the bath 
state are highly symmetric, for instance, if the Hamil- 
tonian is spin isotropic and the bath is prepared initially 
in a completely disordered state. Then we found both nu- 
merically and analytically that the exponent of the scal- 
ing with t acquires an additional factor of 2: d cx r 2< >. 
The same was already observed for the UDD sequence 
applied to pure dephasing in Ref. [29|. 

We emphasize that this result is by no means at odds 
with the proofs of universality [15l - [l7l |. The general 
proofs refer to the worst case for decoherence. They de- 
termine whether a certain operator (Pauli matrix) for the 
qubit occurs or not irrespective of the bath operator to 
which it is multiplied. The underlying idea is that for 
any non-vanishing bath operator there is a bath state 
such that the qubit state is influenced in a non-trivial 
way. Hence decoherence occurs. 

But for certain choices of the bath state even a non- 
vanishing bath operator may have a vanishing effect on 
the quantum bit if its partial bath trace vanishes. Then 
no decoherence is induced by this particular term. This 
is the effect which enhances the performance of the QDD 
sequence for highly symmetric situations. We summarize 
that min{N x , N z } + 1 is a lower bound for £. 

We stress that the found phenomenon is relevant for 
realistic situations. Complete or partial spin symmetry 
in the Hamiltonian is a standard feature. A completely 
disordered bath state is also an excellent starting point 
in the description of baths of nuclear spins. Their mutual 
interaction is so small in energy that even room temper- 
ature suffices to disorder the nuclear spins completely. 

In general, we conclude that the more asymmetric the 
bath Hamiltonian, its coupling to the qubit, and in par- 
ticular its initial state are, the lower the exponent C is of 
the leading non- vanishing power in the total duration r of 
the sequence inducing decoherence. The same is true of 
UDD sequences for pure dephasing. So far, we focussed 
on spin baths which allow for completely disordered, infi- 
nite temperature states. It is an interesting question for 
future research whether a similar phenomenon can occur 
in other baths such as bosonic ones. 

Experimental research is also called for. To our knowl- 
edge, there exist studies on the influence of the initial 
state, see for instance Ref. l32l . but they focus on the ini- 
tial state of the system. Discussions of the influence of 
the initial state of the bath, which was our focus here, arc 
scarce (33j . Moreover, it must be distinguished between 
studies of iterated cycles of sequences with exponential 
decay rates [13] and studies of a single sequence display- 
ing decoherence with a particular power law 0-0, EH, US - 

MM- 

Of course, it is difficult to measure the exponents di- 
rectly. But we suggest to demonstrate experimentally 
that the performance of a QDD or a UDD sequence is 
lowered if the symmetry of cither the Hamiltonian or of 
the initial bath state is lowered. This would already be 
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a smoking gun evidence for the essence of the present 
theoretical finding. 
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Appendix A: Magnus expansion of the evolution 
operator 

In general, it is complicated to find a universal argu- 
ment that explains the scaling of QDD if the symmetry 
of the system is not well defined. The reason why some 
exponents deviate from the analytically predicted value 
depends on the possible vanishing of some terms in the 
Magnus expansion [H], HH of the evolution operator. The 
form of such terms strongly depends on the form of the 
Hamiltonian and of the symmetry of the bath state. 

If no general conclusions can be drawn on the parity of 
these terms and of the density matrix under a rotation 
Pp, one way to proceed is to calculate the cumulants 
of the expansion explicitly and to analyze which terms 
determine the power law of d with r. Here wc provide an 
analysis of the first cumulants of the Magnus expansion 
for N x odd for the data of Tab. IMa). 

For a generic instant t S [0,r] we write the Hamilto- 
nian @ as 

H(t) = H B + ]T o-„/ M (t)V (Al) 

The switching functions [36| f^(t) are the effect of the 
stroboscopic pulsing of the qubit. They are piecewise 
constant functions with values ±1 

/,(*) = (-iy* for te(tj,t? +1 ], (A2a) 

f y (t) = (-l) k for t € (q k ,t* +hk ], (A2b) 

with f y (0) = 1 and f y (r) = —1, and 

Mt) = f z (t)f y (t). (A2c) 

The evolution operator can be written in terms of the 
cumulants as 

oo 

U(t,0) = exp{-ir^iJ ( ™)} (A3) 
n=l 

with rij'"' oc T n . The first and the second cumu- 
lants are defined as tH^ = /J" H(t)dt, tH^> = 

~4? Io dtl Jo' &2 [11(h), H(t 2 )}. From Eqs. JMJ it is 
straightforward to verify that di/ M (t) = for /i = x, 



y or z and for N z — N x — 1. The first cumulant is 
proportional to the bath Hamiltonian 

tH {1) =tH b . (A4) 

For the second cumulant one finds 

2itH^= J2 ^[Hb^] (A5) 

fi—x,y,z 

+ E IfW^avA,,}, (A6) 

fi,u—x 7 y 7 z 

with the integrals 

1% := f dti f 1 dt 2 (/„(i 2 ) - Uih)) (A7a) 
Jo Jo 

and 

I£'":= f T dh f 1 dtaf M {t 2 )Mh). (A7b) 
Jo Jo 

The integrals (jA7a| and ()A7b[) can be easily evaluated, 
we report below only those that are different from zero. 
For N z = N x — 1 one finds 

i! = ^, ii=Y> ir=-ir=~- (as) 

The second cumulant then becomes 

2 2 

2ir #(2) = T —a v [H B ,A y ] + T —a z [H B ,A,] 
t 2 

+ i—a y [A x ,A z } + , (A9) 

where the notation [ , ]+ stands for an anticommutator. 
It is interesting to notice that does not contain any 
terms proportional to <j x . 

If the Hamiltonian is SU(2) invariant and the coupling 
constants Jq and J\ are equal and independent of the 
indexes i and j (J J = aX and J{ = A) the commutators 
in Eq. (|A9|) vanish because the Pauli matrices anticom- 
mute. This is true for a central spin model. 

On the other hand, if A^ = (e.g. for a 

spin chain) the commutators in Eq. (|A9|) are different 
from zero. The anticommutator becomes [Aa;,A z ] + = 

2 iJ2i,j{Ji)xx(Ji)zz( 1 - 5ij)(Tx which implies that it 
vanishes for a spin chain because the qubit is coupled to 
a single site only, i.e., i — j holds. 

As usual we are interested in the difference between 
the evolved density matrix p^\t) = U(t, 0)p 7 ^?7(t, 0)* 
and the initial one. We find 

Ti'B {ti\r) p^} = Tr B {p B [A X ,A Z } + } 

+ ^ c ^y Tr B {[A X1 A Z ] +PB [A X ,A Z } + } (A10) 

where the hermitecity Aj^ = A^ was used. In writing 
Eq. (|A10[) we neglected the contributions coming from 
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the first cumulant because they do not alter the qubit- 
operator content of the norm distance. 

If pb oc 1b the terms with the trace over the bath 
vanish for a SU(2) invariant Hamiltonian (see Sect. IIIip 
while they are finite for an asymmetric model such as the 
one in Eq. © in Sect. El 

In order to understand why the distance norm scales 
with exponents N z + 1 for A^ = 1 (and in general for 
N x odd) in Tab. IPvT ah some knowledge on the third 
cumulant is required. This cumulant is defined (34j 
as 

- f T r* 1 f* 2 

- 6ri? (3) = / dti / dt 2 / dt 3 x 
Jo Jo Jo 

{[H(t 3 ), [H(t 2 ), Hih)}} + [Hfe), [H(t 2 ), H(t 3 )}}} . 

(All) 

The commutators [H(t 2 ), H (tj)] and [H (t 2 ) , H (t 3 )] have 
the same operator content as H but differ in their pref- 
actors and in their time dependence. One can verify that 
[Hb, H^] = if the bath Hamiltonian is SU(2) invariant 
due to the anticommutation of the Pauli matrices. For 
the same reason one finds that 



/ 4 c M,y,f 

/J . IS 



(A12) 



where i ^ j ^ k. Equation (|A12[) only provides the 
operator contents of Eq. (|A11[) . In order to eliminate the 



time dependence one must substitute Eq. (|A12I) into Eq. 
(jAllj) and integrate over ti, t 2 and £3. Each operator 
brings a switching function / M (t) with it, such that the 
integration in (| Al 1|) yields the coefficients 



it^ 1 = r dh r dt 2 r dhf a (h)f (t 2 )f^t 3 ). 

Jo Jo Jo 

(A13) 

The indices a, /3 and 7 can be equal to x, y, and z. We 
have checked numerically that for N x = 1 and N z = 2 
the only non-zero contributions are given by L 



x,z,z j-Z,X,Z 
1 1 3 



and I 3 ' ' corresponding to the qubit operator a x . Thus 
the third cumulant is a term of pure dephasing that 
can be suppressed by means of a Z sequence of 7r pulses. 
From the numerical results we expect that the same ar- 
gument holds in general for higher cumulants and for A^ 
odd. 

We also checked our results for a spin chain. We find 
the same results as for a central spin model in the cases 
of low symmetry, high symmetry and in the case of an 
asymmetric Hamiltonian with pb oc 1b- Discrepancies 
are found for the SU(2) invariant Hamiltonian in combi- 
nation with the product state |^b) ©■ A possible expla- 
nation is provided by the commutators in Eq. (|A9j that 
do not vanish for a spin chain with = . Hence 
the precise topology of the model matters for the case of 
mixed degree of symmetry. 
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